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Abstract 



Generalised observables (POM observables) are necessary for rep- 
resenting all possible measurements on a quantum system. Useful al- 
gebraic operations such as addition and multiplication are defined for 
these observables, recovering many advantages of the more restrictive 
Hermitian operator formalism. Examples include new uncertainty re- 
lations and metrics, and optical phase applications. 

1 Introduction 

The assertion that all observable quantities of a quantum system can be rep- 
resented by Hermitian operators acting on the Hilbert space of the system is 
now well Icnown to be inconsistent with the usual representation of measure- 
ment via interaction between the system and an apparatus. In particular, 
consider the measurement procedure defined by preparing an apparatus in 
some fixed initial state, allowing an interaction between the system and ap- 
paratus, and measuring some Hermitian operator on the apparatus. The 
probability of a given measurement result, a, for initial state IV") of the 
system, is then found to have the general form 



where Aa is a positive operator acting on the system Hilbert space, deter- 
mined by the specific measurement setup (i.e., by the initial apparatus state, 
the interaction Hamiltonian, and the apparatus observable). 

The set of operators {Aa} corresponding to such a measurement proce- 
dure is called a probability operator measure (POM) or a positive operator 
valued measure (POVM), and the corresponding physical observable of the 





system, A, may be referred to as a POM observable ^ [21 El ID- In gen- 
eral Aa is not a projection operator associated with an eigenspace of some 
corresponding Hermitian operator. 

Thus, to consider aZ/ possible measurements on a quantum system (with- 
out reference to specific experimental setups), one must consider all POM 
observables on the Hilbert space of the system. Such observables repre- 
sent a non-trivial extension of Hermitian observables, being necessary to de- 
scribe, for example, optimal measurements for distinguishing between non- 
orthogonal states Eim El ) and for describing optical phase (211110111 
IHl m and optical heterodyne detection ^1 . 

However, what is gained in generality is lost in algebraic simplicity. For 
example, it is not a priori clear how to add and subtract, much less multiply, 
POM observables. It would therefore be desirable to be able to define an 
algebra on the class of POM observables, to regain some of the advantages 
of the operator algebra associated with (the less general) Hermitian observ- 
ables. The aim of this Letter is to demonstrate the existence of such an 
algebra, and applications thereof. 

In the following section POM observables are briefly reviewed, and max- 
imal and non-redundant observables defined. In section 3 the sum, product 
and other binary operations are defined for the relatively simple case of one 
POM observable and one Hermitian observable, and examples of statistical 
deviation and distance given in section 4. In section 5 the algebraic combi- 
nation of two arbitrary POM observables is considered. This general case 
is more difficult, as certain consistency conditions must be observed, but is 
solvable. A generalised uncertainty relation and optical phase examples are 
given in section 6, and conclusions in section 7. 

2 Generalised observables 

The requirement that the probability distribution p{a\'ip) in Eq. is posi- 
tive and normalised for all states ^ implies that the POM {^4^} satisfies 

Aa>0, E^'^ = l- (2) 

a 

In fact, any set of operators {Aa} satisfying these conditions may be realised 
by a measurement procedure of the type discussed in the Introduction, and 
hence Eq. ((SJ fully characterises the class of POM observables 0(21011]. 
The summation is replaced by integration over continuous ranges of mea- 
surement outcomes. 
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From Eq. the expectation value of any function f{A) of a general 
POM observable A is given by 

{f{A)) = J2f{a)p{a\^) = {i;\7{Am, (3) 

a 

where f{A) is defined to be the operator 

JlA):=J2fia)A,. (4) 

a 

Hermitian observables correspond to the special case AaAa' = for a ^ a' 
(and real- valued outcomes). For such observables the associated operator 
j4 = ^ is Hermitian, and satisfies f{A) = f{A) for any function /. Further, 
one has Aa — A^ = Aa J2a'j^a ^a! = 0, i.e., for Hermitian POMs each A^ is a 
projection operator (onto an eigenspace associated with eigenvalue o of A). 

Now, for any positive operator Aa appearing in some POM, a new POM 
can be obtained by replacing Aa by copies of Aa/N (i.e., by ^a,i, • • • , Aa^N-, 
with Aa^i = Aa/N). However, this new POM observable is trivially mea- 
sured, by making a measurement of the original observable and throwing an 
A^-sided die if outcome a is obtained (to distinguish the new possibilities). 
Thus this new POM is physically redundant, only differing from the original 
in a trivial statistical sense that is independent of the actual system. A 
similar redundancy is obtained if the weightings 1/A^ are replaced by any 
set of positive numbers wi, . . . ,wn which sum to unity. 

A POM observable A is therefore called non-redundant if Aa ^ XAf^ for 
a ^ b and any real number A. Note that any POM observable may be 
trivially (and uniquely) reduced to a non-redundant POM observable, by 
adding together any proportional terms in the corresponding POM. Note 
also that non-redundant observables exclude the trivial possibility Aa = 0. 

Finally, a POM observable A is called maximal if Aa = \a){a\ for all 
a (the kets {\a)} are not necessarily orthogonal nor normalised). It may 
be shown that maximal POMs are maximally informative, in the sense that 
any measurement which optimises information gain under a given constraint 
is equivalent to the measurement of some maximal POM. Further, since 
any positive operator Aa may be decomposed into a sum of "back-to-back" 
kets, it follows that every POM observable has a (non-unique) maximal 
extension. Attention will therefore be restricted to non-redundant maximal 
POM observables in what follows. 
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3 Combinations of POM and Hermitian observ- 
ables 

Defining the algebraic combination of a maximal POM observable with any 
Hermitian observable is relatively straightforward, and hence worth exam- 
ining separately from the general case. The key here is the promotion of 
the POM observable to an equivalent Hermitian operator on an extended 
Hilbert space, in such a way that all algebraic relations between existing 
Hermitian operators are preserved. 

In particular, let H denote the Hilbert space of the system, and for a 
maximal POM observable A = {|a)(a|} let i?^ denote the Hilbert space of 
square-integrable functions over the space of outcomes of A. There is then 
an orthonormal basis {|a)} for Ha associated with these outcomes, and a 
natural mapping from H to Ha defined by 

a 

From Eqs. and Q one has 

p{a\^) = \{am^ = \{a\^A)?. 

Thus the POM {[a)(a|} on H has statistics equivalent to the Hermitian 
operator 

A:=^a\a){a\ (6) 

a 

on Ha, where the latter corresponds to the Hermitian POM {[a)(a|}. 

Note that the representation = ("IV') = ("IV'^) is well known for 
particular POM observables (eg, the coherent state and canonical phase 
representations [21 IH])- The associated representation of ^ as a Hermitian 
operator A on Ha provides a simple example of Naimark's extension theorem 
1^ [21 15 for the (admittedly trivial) case of maximal POMs. 

The set of extended states, {IV"^)}, is characterised by the projection 
operator 

E:=Y,{a\a')\a){a'\ = E'^. (7) 

a,a' 

In particular, it may be checked that any normalised ket in Ha is gen- 
erated by some physical state IV'), as per Eq. ©, if and only if 

E\ij) = IV). 
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Thus the unit eigenspace of E is isomorphic to the physical Hilbert space 
H. Note that for a Hermitian observable one has E = 1, and hence in this 
case (and only in this case) Hj^^ is isomorphic to H. 

Further, let X be any Hermitian operator on H. Then there is a natural 
extension of X to a Hermitian operator on Hj^, defined by 

XA = J2{a\X\a)\a){a\. (8) 

a,a' 

It may be checked that (i) the product XY is mapped to X_aY_a, and (ii) 
the state X\ip) is mapped to X^IV'^). Thus this extension preserves all 
algebraic properties of Hermitian operators and states on H. 

Finally, there is an inverse mapping of states and observables from Hj^ 
to H, generated by the projection \a) \a). It follows from Eqs. ©, 
© and Q that \iP_a) — > IV') and X_a X under this mapping, and that 
the Hermitian POM {|a)(a|} corresponding to A is mapped to the original 
POM {|a)(a|}. Moreover, an arbitrary POM {|A;)(A;|} on Hj^ is mapped to 
the POM observable on H, where 

\k):=J2{a\k)\a). (9) 

a 

Note that orthonormality of the basis {|a)} implies 

J2\k){k\ = 5: (a|fc)(fcla')|a)(a'| 

k k,a,a' 

= J2^a\a)\a){a\ 

a,a' 
a 

as required by Eq. ((2)). 

The tools for defining algebraic combinations of A and X (such as the 
sum A + X, the symmetric product {AX + XA)/2, and the "commutator" 

X]), are now laid out. In particular, let g denote any function which 
maps pairs of Hermitian operators to a Hermitian operator. Hence, 

K:=g{A,XA) (10) 

is a Hermitian operator on with an associated Hermitian POM {|/c)(/c|}. 
The corresponding mapping g{A,X), for a general POM observable A and 
Hermitian observable X on H, is then defined to be the POM observable 

g{A,X) ■.= K^{\k){k\}, (11) 
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where |A;) is defined by Eq. @. 

Note from Eqs. © and © that {k\ijj) = {k\'il)X)- Hence the probabihty 
distributions and p{k\il)X) are equivalent, and from Eqs. (fTUl) and ((TT|) 

one therefore has the identity 

{g{AX)) = {^A\g{A,Xj^)\^j^) = Y,{^\a){a\g{A,X_A)\a'){a'\^) (12) 

a, a' 

for expectation values. This is a very useful formula, as it means one can 
calculate the expectation value of g{A, X) without having to explicitly de- 
termine the corresponding POM {k){kW (which would require diagonalising 
the Hermitian operator K = g{A,Xj\,)). 

4 Examples: deviation and distance 

From Eq. (|12j) one may calculate the statistical deviation between A and X, 
for a given state 1-0), as 

{{X-Af) = {^ja\xI + A^-Xj^A-AXj^\^j^) 

= (X^) + (A^) - [{iJA\a){a\Xj,\a'){a'\A\i;j,) + c.c.^ 

a,a' 

= mX-Af\^) + {^p\A^-(Af\^), (13) 

where Eqs. ©, © and (jSJ have been used, and the operators A, A? are 
defined via Eq. Q . This quantity provides a measure of the degree to which 
the POM A may be approximated by a given Hermitian operator X. Note 
that the second term is not equal to the variance of A. One also has the 
alternative formula 

{{X-Af) = 5](V^^|(X^-i)|a)(a|X^-i|V^) 

a 

= Y.Mx-am\ 

a 

which has been used previously (without algebraic justification) in the con- 
text of obtaining exact uncertainty relations between photon number and 
optical phase [12;. 

Further, the "distance" between a POM observable A and a Hermi- 
tian observable X on may be defined via a natural generalisation of the 
Hilbert-Schmidt metric: 

d{A,Xf ■= tr\{A-XAf\_ 

= tx\{X -^f ^ A^ - {A)\ (14) 
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This generalised metric satisfies the triangle inequalities 

d{A, X) + d{A, Y) > d{X, Y) > \diA, X) - d{A, Y) \ 

as an automatic consequence of the Hilbert-Schmidt metric on /7_4, and 
hence d{A, X) may indeed be interpreted as a measure of "distance" . 
From Eq. 1)141) one has the Pythagorean relation 

d{A, Xf = d{A, Af + d(A, Xf 

between A^ X and A. It follows immediately that A represents the Her- 
mitian operator "closest" to A, being separated from A by the minimum 
distance 

mind(AX) = d{A,A) = {iT\^ - (Af]Y/'^ . 

This minimum distance is a useful measure of the inherent "fuzziness" of A, 
vanishing if and only if ^ is a Hermitian observable. 

5 Combining arbitrary POM observables 

To generalise the above results to algebraic combinations of two maximal 
POMs A = {|a)(a|}, B = one must find a Hilbert space if^e which 

contains two suitable orthonormal sets {|a)}, {|6)}. It turns out that, to 
avoid such undesirable properties such as ^ — ^ 7^ 0, the compatibility of 
A and B must explicitly be taken into account. The resulting simultaneous 
mapping of A and B to Hermitian operators on -fT^g corresponds to a novel 
and highly nontrivial Naimark extension, in contrast to the mapping from 
H to H_A defined in section 3. 

First, note that each state \ip) on H will have two possible extensions 
on if^B, given by \'4^jC) and \'4>b) as calculated via Eq. ©, corresponding 
to respective projection operators Eji^ and E'g defined as per Eq. Con- 
sistency requires that these are equal. Calculating {h\il)j{) and (alV";?) via 
Eq. then yields the conditions 

|a)=E(%)l^)' \h)=Y.i^\h)\a) (15) 

6 a 

on the inner product 

Second, since the statistics of A and B on H will be equivalent to those 
of the Hermitian operators A and B on Hji^Q respectively, as calculated via 
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Eq. @, it is physically desirable that any identical statistical components 
of A and A are mapped to identical components of A and B, i.e., 

\a)ia\ = \b){b\ for \a){a\ = \b){b\. (16) 

This condition ensures that compatible aspects of the observables A and B 
are preserved by the extension to H_/(b (and in particular that A = B for 
A = B). 

If conditions (jlSj) and (jlOj) can be satisfied, then any function g which 
maps pairs of Hermitian operators to a Hermitian operator will define a 
Hermitian operator 

k:=g{A,B), 

on Hj[is, with an associated Hermitian POM {|A;)(A;|}. The corresponding 
POM observable g{A, B) on H is then defined to be the POM K = {|A;)(A;|}, 
with 

\k) ■.= Y,{a\k)\a) =Y,{b\k)\b) 

a b 

in analogy to Eq. Q, where the second equality follows immediately from 
Eq. 

It therefore remains to find orthonormal sets {|a)}, {!&)} satisfying Eqs. (|15|1 
and (HH). First, let C := {|a)(a|} n = {\c){c\} denote the set of com- 

mon elements of the POMs corresponding to A and B. C is physically well- 
defined for non-redundant POMs, where measurement of the POM CU{Cq}, 
with 

Co:=l-^Ic)(c|, 

c 

corresponds to the simultaneous measurement of the compatible components 
of A and B. Condition ()16|) implies that these compatible components are 
promoted to simultaneous eigenstates of A and B. 

Decomposing the POMs for A and B as {|c)(c|, |a)(a|} and {|c)(c|, 
respectively, it follows that 

J2\~^){a\=J2\b){b\ = Co. (17) 

a b 

Eqs. (|15|) and (|16j) then imply that the orthonormal set {jc)} is orthogonal 
to each of the orthonormal sets {[a)}, {|^)} on Hj^is, and further that the 
inner product (a|6) satisfies 

\~a)=J2m\b), \b)=Y.{a\b)\~a). 
b 
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These equations are formally solved by the choice 

{a\b) = {d\C,'\b) (18) 

as may easily be checked using Eq. ()17|). However, it must be verified that 
this solution has all the properties required of an inner product between two 
orthonormal (not necessarily complete) sets in some Hilbert space Hj^^q (one 
may then take Hj^^ = H_ab © He, where Hq is the span of {|c)}). 

First, each of \a) and \b) are orthogonal to the zero-eigenspace of Co 
(since Eq. H17|) implies that (a|Co|a) and (6|Co|6) are strictly positive), and 
hence the righthand side of Eq. ()18() is well-defined. Second, noting Cq is 
Hermitian, then {d\b) = {b\d)* as required. Finally, for it to be possible to 
write each \d) as an orthogonal superposition of the \b) and some orthonor- 
mal set {|a^)}, i.e., 

\~a)=Y.m\b) + Y.{x\~a)\x), 
b ^ 

one must have J2b 1(^1^)1^ ^ 1 (and similarly J2a I (^1*^)1^ ^ 

J2 \(b\d)\^ = (a|C7o-i|a) = {a\a) 
b 

1/2 

from Eq. (fTH)) . where |o) := Cq \d). Noting from Eq. (|17() that X]al^)(^l 
is equivalent to the unit operator on the span of {|a)}, one then has 

{a\a) = {a\[J2\a){a']\a) 

a' 

= 5:Ka|a')P>(a|a)2, 

a' 

implying (o|a) < 1 as required. 

6 Examples: uncertainty relations and phase 

The expectation value (and hence the statistics) of any algebraic combina- 
tion of two arbitrary POMs A and B on H may now be calculated via the 
relation 

{g{A,B)) = {ij\g{A,Bm (19) 
analogous to Eq. (fT^ . where 1-0) denotes either of \iPa)i iV'i?)- 
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The examples of statistical deviation and distance in section 4 generalise 
immediately to such pairs of POM observables. As a further example, the 
expectation of the "commutator" of A and B follows as 

{z[A,B]) = i{^|;\[A,B]\^l,) 

= i{ip\[Y,ab(a\b)\a){b\-h.c.]\4)) 

a,b 



i(V'|[E«^(«|Go"'l^)l«)(^l-MIV') 

d,b 

i{^P\AC,'B-BC^'A\i;), 



where the operators 



A:=^a\d){a\, B ■.= Y,h\h){h\ 
a I 

represent the restrictions of A and B to non-identical \a){a\ and \h){h\. The 
Heisenberg uncertainty relation on H_At3 therefore leads immediately to the 
generalised uncertainty relation 

> ]^\{il)\AC^^B - BCq^A\'4))\ (20) 

for arbitrary POM observables A and B on H. 

In the case that A and B have no POM components in common, then 
Co = 1 and the above uncertainty relation reduces to 

A^A0>i|(V'p,F|V)|. 

As an example, for the photon number operator N = X]n'^l^)(^l 
canonical phase POM $ = {|(/')(</'I}, where |(/>) = (27r)-i/2 e*"'^|n) and 
4> E (-vr,7r] [Il|llllEl[71|Hlin],a straightforward calculation yields 



AA^A$ > (l/2)|(V'|[iV,$] 

= (l/2)KV^|^(-ir+"|m)(n| 

m,n 

= (l/2)|l-27rp(7r|V')|, 

in agreement with previous (non-algebraic) methods [HI I12| ITH]. 

Another phase observable of interest arises from ideal heterodyne de- 
tection, where one makes a simultaneous (but noisy) measurement of the 
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quadratures of a single- mode optical field |^ ^1 . This measurement is 
represented by the the so-called "coherent-state" POM Ah = {'^~^\(^){'^\}^ 
where |a) denotes the coherent state corresponding to eigenvalue a = re^^ of 
the photon annihilation operator a. Thus '■= arg^j^ is a corresponding 
phase observable: the heterodyne phase El E] (see also Ref. for a 
different realisation of this observable) . 

Since phase is a periodic variable, the circular deviation 

Sh := 1- |(e**e-^*«)| 

provides a more natural measure of comparison for $ and than does 
(($ _ ^hY) (even better is 1 - |(e*(*~*«))|, but this is more difficult to 
evaluate). This quantity will be close to zero in cases where heterodyne 
phase provides a good approximation to the canonical phase, and close to 
unity in cases where heterodyne phase is a poor estimate of the canonical 
phase. Since 5h is defined by an algebraic function of two POM observables, 
it can be calculated by the above methods. 

In particular, from Eq. (|T!?|) one has {AB) = {'4>\A B\'4j) when A and A 
share no common components. A straightforward calculation gives 

i^* =^|n)(n+l|, 

n 

while a Gaussian integration yields ) 

= ^r(n + 3/2)(n!)-^(n + lY^/'^\n + 

n 

and hence one finds 

6h = 1- Y.\{nm^T{n + 3/2) (n!)-^^ + ^Y^'^- (21) 

n 

The circular deviation between the canonical and heterodyne phases is 
therefore completely determined by the photon number distribution of the 
state. Further, Stirling's formula for the Gamma function may be used to 
obtain the asymptotic formula 

5h ~ (V|(iV + l)"V)/8 

from Eq. H21() . to first order in 1/(A^ + 1). Hence 5h is typically small for 
high energy states, implying that heterodyne detection provides an accurate 
estimate of the canonical phase for such states. 
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7 Conclusions 



It has been shown how to define algebraic operations for arbitrary pairs 
of generahsed quantum observables. As well as being of intrinsic interest, 
this allows many of the advantages of calculations with Hermitian operators 
to be realised for the more general POM observables, as indicated by the 
examples in sections 4 and 6. 

The algebraic combination of any two POM observables A and B with 
any number of Hermitian observables X,Y, Z, . . . may also be defined, with 
Eqs. ((T^ and ((TU]) generahsing to 

{giA, B, X,Y,Z,...)) = B, Xab, Yab, Zab, • • -Mab), 

where \tPab) denotes either of {iPa) = \iPb) and Xj^js denotes either of Xj^ = 
Xb (the second equivalence follows from the first, via E_a = Eb)- Note that 
for the special case where A and B correspond to two Hermitian operators A 
and B respectively, then Hj^is is isomorphic to H, and so g{A, B, X,Y, Z, . . .) 
becomes equivalent to the Hermitian operator g{A, B,X,Y, Z, . . .). 

A number of issues remain for future investigation. First, only binary 
combinations of POM observables have been considered. It is not clear 
whether general combinations of three or more such observables can be con- 
sistently defined, nor even whether, for example, the operation of addition 
is associative. Second, the question of uniqueness has not been examined. It 
is possible there are other solutions satisfying conditions (|15|) and (|16|) . and 
even that one could satisfactorily replace the second of these conditions with 
a weaker (and smoother) requirement along the lines that A is "close" to B 
whenever A is "close" to B. Third, it might be possible to use the existence 
of an algebra for POM observables to provide a "cleaner" formulation of 
QM, not requiring any a priori distinction between POM observables and 
Hermitian observables (which, for example, the usual mapping between clas- 
sical and quantum Hamiltonians requires). For example, the optical phase 
observable can now simply be defined algebraically, in direct analogy to the 
classical formula, as the combination e** := N~^^'^a. 

Applications of the generalised statistical deviation and distance, dis- 
cussed in section 4, to determining the optimal estimate of an observable 
from a given measurement and to finding general "disturbance" and "joint- 
measurement" uncertainty relations, have recently been given (since this 
paper was first prepared) ^3 1^ . 
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